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ABSTRACT 

We have studied the aperiodic variability of the 401 Hz accreting millisecond X-ray pulsar 
SAX J1808.4-3658 using the complete data set collected with the Rossi X-ray Timing Explorer over 
14 years of observation. The source shows a number of exceptional aperiodic timing phenomena that 
are observed against a backdrop of timing properties that show consistent trends in all five observed 
outbursts and closely resemble those of other atoll sources. We performed a detailed study of the 
enigmatic ~410 Hz QPO, which has only been observed in SAX J1808.4-3658. We find that it ap¬ 
pears only when the upper kHz QPO frequency is less than the 401 Hz spin frequency. The difference 
between the ^410 Hz QPO frequency and the spin frequency follows a similar frequency correlation 
as the low frequency power spectral components, suggesting that the ~410 Hz QPO is a retrograde 
beat against the spin frequency of a rotational phenomenon in the 9 Hz range. Comparing this 9 Hz 
beat feature with the Low-Frequency QPO in SAX J1808.4-3658 and other neutron star sources, we 
conclude that these two might be part of the same basic phenomenon. We suggest that they might be 
caused by retrograde precession due to a combination of relativistic, classical and magnetic torques. 
Additionally we present two new measurements of the lower kHz QPO, allowing us, for the first time, 
to measure the frequency evolution of the twin kHz QPOs in this source. The twin kHz QPOs are seen 
to move together over 150 Hz, maintaining a centroid frequency separation of (0.446 ± 0.009)r'spin- 
Subject headings: pulsars: general - stars: neutron - X-rays: binaries - individual (SAX J1808.4- 
3658) 


1. INTRODUCTION 

The transient X-ray source SAX J1808.4-3658 
(SAX J1808) was d iscovered with the Be ypoSax satellite 
in September 1996 (in’t Zand et al.|1998 ). It remained in 
quiescence for 1.6 yr alter the discovery, until the Rossi 
X-ray Timing Explorer (RXTE) detected rene wed X-ray 
activity wit h the Proportional Counter Array (Marshall 


et al. I998|). The initial RXTE observations revealed 
coherent pulsations at 401 Hz, making SAX J1808 the 
first known accreting millisecond X-ray pulsar (AMXP) 


( Wijnands fc van der Klis||1998a ). 

I'rom extrapolating the recurrence time between the 
1996 and 1998 outbursts, SAX J1808 was expected to 
show a new cycle of activity around November 1999. 
Unfortunately, Solar constraints prevented observations 
until late January 2000, at which time SAX J1808 was 
observed to be in a low luminosity state at r oughly a 
tenth of the typ ical flux of the 1998 outburst ( 

Klis et al. 2000). It is generally assumed that tlie mam 
outburst occ urred when the Solar constraints prevented 
observation (Wifnands et al.||2001 ). The low luminosit y 
state persisted for ~ 114 days (Wijnands et al. 2001), 
showing dramatic changes in luniinosity. Un a timescale 
of ~5 days the luminosity would drop below a detec¬ 
tion limit of ^ 10^^ erg s“^ and then increase again 
to ~ 10^® erg s“^. During these episodes of increased 
emission the 401 Hz pulsations were detected al ong with 
a sporadic stron g modulation at about 1 Hz (van der 
Klis et al.||2000D. |Campana et al.| (|2008|) proposed that 
the 5-day episodic variations of luminosity w ere caused 
by an intermitte ntly active propeller effect (Illarionov 
& Sunyaev 1975), periodically inhibiting accretion onto 
the neutron star. The 1 Hz flaring remained unique to 
SAX J1808 until a similar feature was discovered in the 


AMXP NGC 6440 X-2 (jPatruno et al.||2010t |Patruno & 
D’Angelo|[2M^ . 


The fourth outburst of SAX J1808 was detected in Oc¬ 
tober 2002, some 2.8 yr after the 2000 outburst, trig¬ 
gering an extensive monitoring campaign with RXTE. 
Twin kHz quasi-periodic osc illations (QPOs) were de¬ 
tected early in the outburst ([Wijnands et al. |2003|), for 
the first time in a source for which the spin frequency was 
precisely known. An additional QPO was unexpectedly 
detected at 410 Hz, which so far is unique to SAX J1808 
and whose nature remains a mystery. After the main 
outburst, which lasted several weeks, the source entered 
a low luminosity state showing the same ~5 day epis odes 
of increased emission and 1 Hz flaring seen in 2000 (Wi- 
jnandspOOd ). 


Since its 2002 outburst, SAX J1808 has been observed 
in outburst another three times; in 2005, 2008 and 2011. 
For both the 2005 and 2008 outbursts a low-luminosity 
prolonged outburst tail was observed, again showing lu¬ 
minosity variations on a ~ 5 day timescale. During these 
episodes the 1 Hz flaring was again observed in 2005 but 
not in 2008 (|Patruno et al. 2C)09|). Observations of the 
2011 outburst were cut short by Solar constraints, such 
that only the onset of the low-luminosity outburst tail 
was seen, where ag ain, no 1 Hz flaring was detected (]Pa-| 


truno et al. 2012). The 2008 and 2011 outbursts were 
later found to show a 1-5 Hz flaring which shares many 
of the characteristics of the 1 Hz flaring, but occurred at 
luminosities an order of magnitude higher, near the peak 
of the outburst (Bult & van der Klis 2014). 


SAX J1808 showed several thermonuclear X-ray bursts 
in each outburst. Using the X-ray bursts of the 1998 , 
2002 and 2005 outbursts Galloway & Gumming (2006) 
derived a distance estimate of 3.5 ± 0.1 kpc. 

The literature contains a large body of work on the co- 
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herent pulsations of SAX J1808 (Wijnands & van der Klis 

I998al fPoutanen & Gierlihski||20U3| Burderi et al. 2006 

Hartman et al. 1120081 Leahy et al. 

120081 Hartman et al. 

2009 Ibragimov & Poutanen|2009 

Patruno et al.||20I2p. 


on the other hand, remains incomplete, with all work 


and the corresponding colors and intensities were consis¬ 
tent with being the same. Details of the intervals used are 
given in Table Finally we estimated the background 
count rate using the ftool pcabackest and used it 
to renormalize the power spectra to units of source frac- 


jnands & van der Klis 

(1998b 

) studied the broadband Lorentzian prohles (see e.g. 

Nowak 

2000 j 


power spectrum of the 1998 outburst and |Menna et al. 


(20031 considered the coupling between the broad band 


noise and the coherent pulsatio ns. The 2002 outburs t 
power spectrum was studied by Wijnands et al. (2003), 
but only the high frequency region, in the context of the 
kHz QPOs. A more gener al analysis of the 1998 an d 
2002 outbursts was done by van Straaten et al. (2005), 
who compared the timing behavior of AMXFs with non- 
pulsating low-mass X-ray binaries. The stochastic vari¬ 
ability of the 2005, 2008 and 2011 outbursts remains un¬ 
addressed. 

With the demise of RXTE, its 14 year monitoring cam¬ 
paign of SAX J1808 has now been completed. We use this 
remarkable dataset to provide a complete overview of the 
stochastic time variability of SAX J1808. In this work 
we analyze all outbursts in a consistent approach, but do 
not consider the low lum inosity outburs t tail, which has 
been studied in-depth by Patruno et al. (2009). We pro¬ 
vide the first in-depth analysis of the broadband power 
spectra of the 2005, 2008 and 2011 outbursts and char¬ 
acterize the rich variability phenomenology observed in 
this accreting millisecond X-ray pulsar. 

2. DATA REDUCTION 

We used all RXTE pointed observations of 
SAX J1808’s outbursts (see Table for ObsIDs), 
selecting only observations with stable pointing, source 
elevation above 10° and pointing offset less than 0.02°. 

Using the Standard-2 data we created Crab normal¬ 
ized 2-16 keV light curves. We also calculated the Crab 
normalized soft color as the count rate in the 3.5-6.0 keV 
band divided by the count rate in the 2.0-3.5 keV band. 
Similarly, the Crab normalized hard color is calculated 
by dividing the 9.7-16 keV band cou nt rate by the count 
rate i n the 6.0-9.7 keV band (see e.g. van Straaten et al. 
2003| for the detailed procedure). 

The timing analysis was done using high time resolu¬ 
tion (122 fj,s) Event and GoodXenon data, selecting all 
events in the 2-20 keV energy band. Data obtained in 
GoodXenon mode were binned to 122 fis prior to fur¬ 
ther analysis. Fourier transforms were calculated using 
256 s segments at full time resolution, yielding power 
spectra with a frequency resolution of ~0.004 Hz and 
a Nyquist frequency of 4096 Hz. No background sub¬ 
traction or dead-time correction was applied prior to the 
Fourier transform. The power spectra we re calculated 
using the standard Leahy nor malization (Leahy et al . 
1983). Following the method of Klein-Wolt et al. 1 20U4| 
we inspected the high frequency region (> 1500 Hz 
for anomalous features. Since none were found we sub¬ 
tracted a Poisson noise p ower spectrum calc ulated using 
the analytical formula of Zhang et al. (1995) and shifted 
to fit the > 1500 Hz power. All power spectra were aver¬ 
aged per ObsID. If necessary, multiple consecutive Ob¬ 
sIDs were averaged into intervals to improve statistics, 
but only if the power spectra showed similar morphology 


h'ourier frequency v and determined by three parameters: 
the characteristic frequency i/max = i/r'o + (FWHM/2)2, 
the quality factor Q = r'o/FWHM and its fractional rms 
amplitude squared = P = L{v)dv. Here vq is 
the Lorentzian centroid frequency and FWHM the full- 
width-at-half-maximum. A feature is called a QPO if 
Q > 2. Features with Q < 2 are called noise compo¬ 
nents. We consider a component to be significant if the 
ratio of the integrated power to its single trial negative 
error is greater than three; P/ap > 3. Boundary cases 
with P ^ iap for which we deviate from this rule are 
explicitly discussed in Section]^ 

In the 2008 and 2011 outbursts we previousl y found 


a broad noise com ponent centered at 1-5 Hz (Bult & 


van der Klis 

2014 

), which we called high luminosity flar- 

mg. 'i’his t 

taring component cannot be litted with a 


Lorentzian profile, but instead is well described by a 
Schechter function, which is a powe r-law with an expo- 


Klis 

1989 

Dotani et al. 

1989 


function with Lorentzian components, and by analogy to 
the definitions used for Lorentzian components, we define 
the Schechter “centroid” frequency, vq = —az/cut: as the 
frequency of maximum power in S{v). Alternatively, we 
can define a characteristic frequency, zzmax = (1 ~ Q;)iZcut, 
as the frequency of the maximum in uS{v). 

3. SOURCE STATE IDENTIFICATION 

Atoll sources are named after the shape they trace out 
in the color-color diagram (Hasinger fc van der Klis|1989 


van der Klis 2006). At high luminosities they exhibit a 
banana shaped pattern, which is divided into three re¬ 
gions corresponding to three source states {upper, lower 
and lower-left banana). At lower luminosities the hard 
color increases and the observations cluster in the island 
state. In the hardest, lowest luminosity states a change 
in variability distinguishes another source state; the ex¬ 
treme island state. Depending on the source state, atoll 
sources show a different number of power spectral fea¬ 
tures and have a different morphology (see van der Klis 
2006 for a detailed overview). 


5AX J1808 is a low luminosity atoll source (van 
Straaten et al. 2005) and has shown only three source 
states: the extreme island state, the island state and the 
lower-left banana. We discuss the power spectral shapes 
relevant to these source states. 

In the extreme island state three to four broad 
Lorentzians are needed to fit the power spectrum. The 
components with the lowest and second lowest character¬ 
istic frequencies are res pectively called the break (6) and 


hump {h) Lorentzians (Belloni et al. 2002|). The highest 
frequency component is usually found at 11)0-200 Hz and 
to represent the widely held idea that 
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Table 1 
Interval Listing 


Interval 

State 

Start MJD 

ObsIDs 


1998 Outburst 

1 

IS 

50914.9 

A-01-03S 

2 

EIS 

50919.8 

A-03-00 

3 

EIS 

50920.1 

A-04-00 

4 

EIS 

50921.3 

A-05-00 

5 

EIS 

50921.5 

A-06-01 

6 

EIS 

50921.7 

A-06-000 

7 

EIS 

50921.9 

A-06-00 

8 

EIS 

50923.9 

A-07-00 

9 

EIS 

50926.8 

A-08-00 

10 

EIS 

50927.7 

A-09-01, A-09-02 

11 

EIS 

50928.6 

A-09-03, A-09-04 

12 

EIS 

50929.8 

A-09-00 

13 

EIS 

50930.6 

A-10-02, A-10-01, A-10-00 

2002 Outburst 

1 

IS 

52562.1 

B-03-04-00 

2 

LLB 

52562.3 

B-01-01-000 

3 

LLB 

52563.1 

B-01-01-03, B-01-01-04 

4 

LLB 

52563.4 

B-01-01-01 

5 

LLB 

52564.3 

B-01-01-020, B-01-01-02 

6 

IS 

52565.0 

B-01-02-01, B-01-02-000, B-01-02-00 

7 

IS 

52566.0 

B-01-02-02, B-01-02-03, B-01-02-04 

8 

IS 

52567.1 

B-01-02-05, B-01-02-06, B-01-02-07 

9 

IS 

52568.1 

B-01-02-10, B-01-02-08 

10 

IS 

52569.0 

B-01-02-20 

11 

IS 

52569.2 

B-01-02-09 

12 

IS 

52569.4 

B-01-02-19, B-01-02-23 

13 

IS 

52570.0 

B-01-02-21, B-01-02-11 

14 

IS 

52570.4 

B-01-02-18, B-01-02-12, B-01-02-22 

15 

IS 

52571.0 

B-01-02-13, B-01-02-15, B-01-02-14, B-01-02-16 

16 

IS 

52572.0 

B-01-02-17, B-01-03-00, B-01-03-04, B-01-03-05 

17 

IS 

52572.2 

B-01-03-000 

18 

IS 

52572.9 

B-01-03-06, B-01-03-07, B-01-03-010, B-01-03-01 

19 

IS 

52573.4 

B-01-03-08, B-01-03-09, B-01-03-10, B-01-03-11, B-01-03-12, B-01-03-13 

20 

IS 

52574.1 

B-01-03-02 

21 

IS 

52574.8 

B-01-03-14, B-01-03-03, B-03-05-00, B-02-01-04 

22 

IS 

52577.1 

B-02-01-000, B-02-01-00 

2005 Outburst 

1 

IS 

53523.0 

C-01-01 

2 

IS 

53523.5 

C-01-02, C-01-03, C-01-04, C-01-05, C-02-00 

3 

IS 

53524.9 

C-02-01 

4 

IS 

53525.8 

C-02-02 

5 

IS 

53526.8 

C-02-030 

6 

IS 

53527.0 

C-02-03 

7 

LLB 

53527.8 

C-02-04, C-02-05, C-02-06 

8 

IS 

53530.8 

C-02-08 

9 

IS 

53531.3 

C-03-03, C-03-02, C-03-00, C-03-01 

10 

IS 

53533.6 

C-03-05, C-03-04, C-03-07, C-03-08 

11 

IS 

53535.7 

C-03-06, C-03-09, C-03-11, C-03-10 

2008 Outburst 

1 

IS 

54731.9 

D-01-01, D-01-00, D-01-02 

2 

IS 

54732.7 

D-01-080, D-01-08, D-01-03 

3 

IS 

54733.7 

D-01-04, D-01-07, D-01-06, D-01-05, D-01-10, D-02-00, D-02-01 

4 

IS 

54736.5 

D-02-05, D-02-06, D-02-03, D-02-04 

5 

IS/F 

54739.4 

D-02-07, D-02-09, D-02-02, D-02-08 

6 

IS 

54742.5 

D-03-00, D-03-01, D-03-08, D-03-02 

7 

IS 

54745.5 

D-03-03, D-03-10, D-03-05, D-03-04 

2011 Outburst 

1 

IS/F 

55869.9 

E-01-01, E-01-00 

2 

IS/F 

55871.1 

E-01-02 

3 

IS/F 

55872.0 

E-01-03, E-01-04 

4 

LLB/F 

55872.9 

E-01-05, E-01-06 

5 

LLB 

55874.0 

E-01-070, E-01-07 

6 

IS 

55875.2 

E-01-08, E-01-09, E-02-00 

7 

IS 

55877.0 

E-02-01, E-02-02 

8 

IS 

55878.0 

E-02-03, E-02-04 

9 

EIS 

55879.0 

E-02-05, E-02-06, E-02-07, E-02-08 

10 

EIS 

55881.3 

E-02-09, E-02-10, E-03-00, E-03-01, E-03-04 


Note. — RXTE ObsIDs grouped together in intervals for all outbursts. ObsIDs are chronologically ordered from left to right and top to 
bottom. Source states are abbreviated as: EIS - extreme island state; IS - island state; LLB - lower-left banana, F - low luminosity flaring present. 
A = ‘30411-01’; B = ‘70080’, C = ‘91056-01’, D = ‘93027-01’, E = ‘96027-01’. 
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this component evolves into the upper kHz QPO found in 


the is land state and the lower-left banana (Psaltis et al. 
1999|). An additional component is usually present be- 


tween the upper kHz and hump Lorentzian, which is pro¬ 
posed to evolve into the lower kHz QPO (iPsaltisetaL 
1999|, but that identification is debated (Ivan btraaten 


et ai. |2003 ). In this work we follow van Straaten et al. 
(|2UU5| and label this feature as a separate “£ow” compo- 
nent. The relation between the (.aw compone nt and the 
lower kHz QPO is further discussed in section [5.1.2| 

In the island state three or four Lorentzian prohles are 
needed. These components tend to be less broad, while 
covering a higher frequency range with respect to the 
Lorentzians seen in the extreme island state. The lowest 
frequency noise components are again called the break 
and hump Lorentzians. The highest frequency compo- 
nent is now a QPO, and identified as the upper kHz QPO 
(Wijnands et al. 2003). Often an additional Lorentzian 
is needed at t'max ~ 100 — 200 Hz, which is called the 


below 0.5. Occasionally it becomes more peaked, but it 
always has Q < 2. A characteristic property of the hHz 
component is that it does not show a systematic trend 
in frequency. It stays roughly constant in the hectohertz 
range. A (ow component is not usually observed in the 
island state, but given their similar frequencies, the (ow 
and hHz components may be difficult to disentangle. We 
simply always label the noise component below the upper 
kHz QPO as the hHz c omponent and discuss this issue 
further in section [5.1.31 

In the lower-left banana the power spectra show a com¬ 
plex morphology and can require up to seven different 
Lorentzians. In contrast to the two island states, the 
lower-left banana may show two breaks in the broad band 
noise, “5” and “&2”, with 62 being the lower frequency 
component. As in the island state, higher frequency fea¬ 
tures are called the hump and hectohertz components. 
Finally, at the highest frequencies a lower {() and upper 
(u) kHz QPO is seen. If both kHz QPOs are detected 
in SAX J1808, the upper kHz QPO is always the most 
prominent feature, so if only one QPO is present, it is 
identified as the upper kH z QPO ( Wijnands et al.|[2003 


van Straaten et al. 2005|). Because the lower kHz QPO 
always appear as a very sharp, weak QPO it can be eas¬ 
ily distinguished from the much broader hHz component 
that may be present at a lower frequency. 

In some atoll sources the hump feature occasion¬ 
ally shows two discernible peaks, or a complex mor¬ 
phology requiring two Lorentzian profiles for an ac¬ 
ceptable fit. In such cases the narrow Lorentzian is 
called the Low-Frequency (LF) component, while the 
broader Lorentzian is called the hump (6) component. 
In SAX J1808 the two components are sometimes both 
narrow QPOs, in which case we call the lower frequency 
component LF and the higher frequency component the 
hump. 

Following 


fBelloni 

profile” 


et 


al. 

"E” 


(20021 we refer to each 
where subscript i = 


Lorentzian pronie as 
{62, 6, LF, 6, £ow, hHz, It} represents the name of the 
specific component (e.g. Lh refers to the break compo¬ 
nent). We may further use Vi to refer to characteristic 
frequency and tzg.i for the centroid frequency. Similarly, 


hectohertz (hHz) component ( 

Ford & van der Klis 

1998 1,1 

van Straaten et al.|2002 

). 'This component is usual 

ly ob- o 

served as a broad noise component with a quality factor “ 



10 15 

Time (days) 

Figure 1. Evolution of the outbursts of SAX J1808. The top 
panel gives the times of type I X-ray bursts. The second panel gives 
the 2—16 keV light curves, while the third and fourth panel give the 
evolution of the soft and hard colors, respectively. Color measure¬ 
ments with large uncertainties have been omitted for clarity. Each 
light curve has been shifted in time by: 2011 MJD 55867.03; 2008 
MJD 54731.20; 2005 MJD 53521.79; 2002 MJD 52556.19; and 1998 
MJD 50906.03, such that their color minima approximately match 
on day 7.0. Each point in this figure corresponds to an ObsID in 
TabletH 

the quality factor and fractional rms are referred to as 
Qj and r^. 

4. RESULTS 

We briefly discuss the timing and color evolution of 
SAX J1808 for each individual outburst. We make use 
of Figure which shows the 2-16 keV light curves, and 
the soft and hard color evolution for all outbursts. Ad¬ 
ditionally, we refer to Figure which shows a sample of 
representative power spectra Tor each outburst, and Ta¬ 
bles and which list the full set of fit parameters 

obtains for SAX J1808. 

4.1. 1998 

The first RXTE detection of SAX J1808 was with 
a scan observation on MJD 50915 (day 9 of Figure 
Marshall et al. 1998), during which the source flux was 
59 mCrab. At this time the source was in the island state 
showing a power spectrum that can be fitted with three 
Lorentzians (if,, Lh and Lu). 

The first pointed RXTE observation was five days later 
on MJD 50920, at which time the flux had decayed to 
^43 mCrab. The soft color was decreasing, while the 
hard color increased. The power spectrum shifted to 
lower frequencies and requires an additional L^ow com¬ 
ponent, clearly showing that SAX J1808 transitioned to 
the extreme island state. 
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Figure 2. Selected power spectra of each outburst of SAX J1808. Every row represents a separate outburst as indicated. Each power 
spectrum has been marked with the number of the interval in its respective outburst and its associated source state (see Table^for details). 
Parameters of the fitted curves are given in Tables [2j|^ 
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Table 2 

Characteristic Frequency; i^max (Hz) 


Interval 

Lb2 

Lb 

Lh ±LF 

-^£ow 

-^hHz 

Li 

Lu 


1998 Outburst 

1 


1.83 ±0.08 

13.3 ±0.5 




236 ± 14 

2 


0.408 ±0.013 

2.89 ±0.08 

26 ±2 



155 ± 7 

3 


0.381 ± 0.018 

2.79 ±0.10 

24 ±4 



167 ± 13 

4 


0.363 ±0.012 

2.62 ±0.07 

21.9 ± 1.8 



150 ±6 

5 


0.38 ± 0.03 

2.39 ±0.16 

22 ±3 



162 ± 20 

6 


0.347 ± 0.013 

2.57 ±0.09 

28 ±3 



177 ± 11 

7 


0.359 ±0.014 

2.63 ±0.09 

21 ± 3 



164 ± 11 

8 


0.37 ± 0.02 

2.80 ±0.17 

27 ± 6 



159 ± 16 

9 


0.81 ± 0.04 

5.5 ±0.3 

54 ±8 



213 ± 14 

10 


1.31 ± 0.07 

7.6 ±0.4 

52 ±8 



199 ± 19 

11 


2.01 ± 0.08 

13.1 ±0.8 




170 ± I4I 

12 


1.03 ±0.03 

6.5 ±0.4 

32 ± 10 



193 ± 34 

13 


0.93 ± 0.03 


37 ± 4 




2002 Outburst 

1 


10.4 ±0.5 

43.9 ±2.1 


179 ± 26 


567 ± 25 

2 


9.9 ±0.2 

57.8 ±2.2 37.4 ±1.1 


212 ± 34 

435 ± 12 

599 ±5 

3 

7.9 ± 1.5 

15.4 ±0.4 

78.2 ±2.0 47.2 ±1.8 


257 ± 35 


698 ±6 

4 

10.0 ± 1.0 

15.9 ±0.2 

75.7 ± 1.0 47.4 ±0.7 


355 ± 43 

504 ± 3 

686 ±5 

5 

2.6 ±0.5 

7.6 ±0.2 

39.3 ±1.1 28.1 ±1.1 


101 ± 13 


497 ± 8 

6 


3.14 ±0.05 

18.01 ±0.15 


122 ±9 


352 ± 5 

7 


3.08 ± 0.05 

18.53 ±0.19 


122 ± 11 


339 ±5 

8 


3.00 ± 0.05 

18.40 ±0.19 


127 ± 12 


332 ± 6 

9 


2.26 ±0.05 

13.6 ± 0.3 


112 ± 12 


294 ±7 

10 


4.6 ±0.2 

24.4 ±1.0 


98 ±48 


389 ± 21 

11 


5.00 ±0.12 

25.4 ±0.3 


113 ± 14 


393 ±5 

12 


4.9 ±0.2 

25.1 ±0.6 


100 ± 25 


381 ± 10 

13 


4.26 ±0.11 

23.0 ± 0.3 


120 ± 19 


370 ±7 

14 


3.41 ± 0.17 

19.3 ±0.5 


89± 28 


332 ± 11 

15 


3.71 ± 0.10 

20.2 ±0.3 


109 ± 10 


358 ± 10 

16 


3.53 ±0.19 

23.2 ±0.7 


88 ± 7 


336 ± 17 

17 


3.41 ± 0.10 

20.9 ± 0.4 


127 ± 22 


326 ± 10 

18 


3.15 ±0.07 

19.2 ±0.3 


139 ± 27 


340 ± 11 

19 


3.49 ±0.14 

20.8 ±0.5 


187 ± 48 


358 ± 14 

20 


3.73 ±0.12 

21.2 ±0.4 


139 ± 18 


332 ± 11 

21 


3.99 ±0.15 

23.1 ±0.5 


173 ± 51 


357 ± 16 

22 


4.8 ±0.4 

28.7 ±1.0 




413 ± 143 

2005 Outburst 

1 


2.6 ±0.2 

23.7 ±1.9 





2 


3.15 ±0.06 

18.9 ±0.2 


122 ± 21 


304 ± 10 

3 


5.51 ± 0.12 

28.1 ±0.4 


90 ± 10 


399 ±7 

4 


7.50 ±0.14 

38.5 ± 0.6 


155 ± 35 


526 ±8 

5 


7.63 ±0.12 

38.4 ± 0.6 


121 ± 21 


532 ± 5 

6 


6.67 ± 0.17 

34.2 ±0.7 


161 ± 45 


494 ± 14 

7 

8.9 ± 1.0 

15.4 ±0.3 

78 ± 4 47.9 ± 0.6 


438 ± 56^ 


685 ± 10 

8 


7.24 ± 0.15 

36.6 ±0.5 


136 ± 12 


498 ± 11 

9 

3.3 ±0.7 

7.9 ±0.3 

36.5 ±0.4 


131 ± 19 


488 ±9 

10 


6.70 ±0.16 

34.4 ±0.4 


127 ± 27 


481 ± 19 

11 


5.07 ± 0.18 

28.3 ± 0.6 


138 ± 30 


393 ± 25 

2008 Outburst 

1 


2.85 ±0.10 

18.7 ±0.8 


135 ± 17 


350 ± 36 

2 


3.66 ±0.17 

19.5 ± 0.4 


118 ±23 


372 ± 14 

3 


5.35 ±0.15 

26.5 ±0.5 


135 ± 29 


415 ± 5 

4 


9.2 ±0.5 

49 ± 2 




524 ± 33 

5 


_ 2 

37 ± 3 




606 ± 44 

6 


8.7 ± 1.3 

54 ±10 





7 


1.72 ±0.16 

23 ± 4 





2011 Outburst 

1 


_ 2 

49.1 ±1.4 


291± 203 


657 ± 3 

2 


_ 2 

38 ± 16 




527 ± 226 

3 


_ 2 

52 ±8 52.3 ±0.7 


343 ± 22 


719 ±4 

4 


_ 2 

59 ±14 51.6 ±1.5 


289 ± 43 


712 ± 5 

5 

4.4 ± 1.1 

15.3 ±0.8 

78 ±4 48.8 ±1.7 



567 ± 9 

737 ± 10 

6 


6.8 ±0.4 

30.3 ±1.0 


153 ± 71 


478 ± 25 

7 


3.7 ± 0.3 

21.1 ±0.9 




245 ± 41 

8 


1.90 ±0.13 

22 ± 3 





9 


0.55 ± 0.04 

4.6 ±1.2 




113 ±36 

10 


0.50 ±0.04 

10 ± 7 






^ Uncertain identification. 

^ Schechter function fit to the high luminosity flaring, see Table 
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Table 3 

Quality Factor; Q 


Interval 

Lb2 

Lb 

Lh 

-^LF 

Lhiii, 

Li 

Lu 


1998 Outburst 

1 


0.37 ± 0.04 

0.29 ±0.07 




0.43 ± 0.10 

2 


0.21 ± 0.03 

0.34 ± 0.06 

0.11 ±0.14 



0.35 ±0.08 

3 


0.19 ±0.04 

0.41 ± 0.08 

0 ( fixed ) 



0.24 ± 0.10 

4 


0.18 ±0.03 

0.41 ± 0.06 

0 ( fixed ) 



0.38 ±0.07 

5 


0.25 ±0.06 

0.43 ±0.12 

0 ( fixed ) 



0.32 ±0.16 

6 


0.22 ±0.03 

0.23 ±0.06 

0 ( fixed ) 



0.30 ±0.10 

7 


0.23 ±0.03 

0.35 ± 0.08 

0 ( fixed ) 



0.21 ± 0.11 

8 


0.28 ±0.05 

0.24 ±0.12 

0 ( fixed ) 



0.41 ± 0.19 

9 


0.20 ±0.03 

0.28 ±0.13 

0.19 ±0.16 



0.8 ±0.2 

10 


0.16 ±0.03 

0.5 ±0.2 

0.4 ±0.2 



0.5 ±0.2 

11 


0.12 ±0.02 

0.57 ± 0.19 




0 ( fixed 

12 


0.15 ±0.02 

1.1 ± 0.5 

0 ( fixed ) 



0.4 ± 0.3 

13 


0.02 ±0.03 


0 ( fixed ) 




2002 Outburst 

1 


0.35 ±0.04 

0.89 ±0.18 


1.2 ±0.7 


0 ( fixed ) 

2 


0.35 ±0.02 

1.9 ±0.7 

2.4 ±0.6 

0.4 ±0.3 

12 ± 113 

7.7 ± 1.7 

3 

0.35 ±0.06 

1.3 ±0.3 

2.8 ±0.8 

2.0 ±0.5 

0.6 ±0.3 


8.2 ± 1.9 

4 

0.28 ± 0.02 

1.3 ±0.2 

3.4 ±0.6 

2.3 ±0.3 

0 ( fixed ) 

27 ( fixed ) 

7.0 ± 1.2 

5 

0.36 ±0.07 

0.58 ±0.10 

1.9 ±0.4 

1.9 ±0.3 

0 ( fixed ) 


3.0 ±0.3 

6 


0.236 ±0.012 

0.85 ± 0.04 


0.19 ±0.08 


2.3 ±0.3 

7 


0.260 ±0.014 

0.83 ± 0.04 


0.20 ±0.11 


2.3 ±0.4 

8 


0.279 ±0.013 

0.86 ± 0.05 


0.24 ±0.12 


2.4 ± 0.5 

9 


0.217 ± 0.016 

0.65 ± 0.06 


0.14 ±0.11 


1.9 ±0.5 

10 


0.33 ±0.04 

1.1 ± 0.3 


0 ( fixed ) 


2.8 ± 1.8 

11 


0.267 ± 0.018 

1.18 ±0.10 


0.15 ±0.19 


4.4 ± 0.6 

12 


0.30 ±0.03 

1.14 ±0.17 


0 ( fixed ) 


5±3 

13 


0.265 ±0.019 

0.93 ±0.07 


0.3 ±0.2 


2.8 ±0.8 

14 


0.27 ± 0.04 

0.98 ±0.16 


0 ( fixed ) 


4± 2 

15 


0.240 ±0.018 

0.99 ± 0.08 


0.32 ±0.16 


2.5 ±0.5 

16 


0.26 ±0.04 

0.72 ±0.10 


1.4 ±0.6 


1.5 ±0.4 

17 


0.255 ±0.021 

0.88 ± 0.09 


0.3 ±0.2 


2.6 ±0.9 

18 


0.245 ±0.017 

0.91 ± 0.07 


0.18 ±0.16 


2.0 ±0.7 

19 


0.25 ±0.03 

0.85 ±0.10 


0 ( fixed ) 


3.7 ± 1.9 

20 


0.24 ± 0.02 

1.01 ± 0.09 


0.34 ±0.19 


3.0 ± 1.3 

21 


0.24 ± 0.03 

1.04 ±0.11 


0 ( fixed ) 


3.2 ± 1.6 

22 


0.21 ± 0.06 

1.4 ±0.2 




0 ( fixed ) 

2005 Outburst 

1 


0.49 ±0.12 

0.7 ± 0.2 





2 


0.315 ±0.017 

1.08 ±0.07 


0 ( fixed ) 


3.4 ± 1.7 

3 


0.324 ± 0.019 

1.28 ±0.10 


0.40 ±0.16 


4.5 ±0.9 

4 


0.308 ±0.016 

1.31 ± 0.11 


0 ( fixed ) 


4.4 ± 1.5 

5 


0.319 ±0.014 

1.37 ± 0.12 


0 ( fixed ) 


5.8 ± 1.9 

6 


0.41 ± 0.03 

1.6 ±0.2 


0 ( fixed ) 


9±6 

7 

0.26 ±0.02 

1.3 ±0.2 

0.73 ±0.15 

4.0 ±1.0 

0.8 ± 0.51 


7±6 

8 


0.293 ±0.018 

1.63 ±0.12 


0.53 ±0.18 


3.8 ± 1.0 

9 

0.34 ± 0.04 

0.69 ±0.11 

1.30 ±0.08 


0.24 ±0.15 


2.7 ± 0.5 

10 


0.253 ±0.019 

1.56 ±0.12 


0.3 ±0.3 


2.3 ±0.9 

11 


0.18 ±0.02 

0.95 ±0.11 


0 ( fixed ) 


2.1 ± 0.8 

2008 Outburst 

1 


0.35 ±0.04 

0.79 ± 0.09 


1.0 ±0.4 


2.4 ± 1.0 

2 


0.20 ±0.04 

1.45 ±0.18 


0 ( fixed ) 


4± 4 

3 


0.28 ±0.02 

1.20 ±0.11 


0 ( fixed ) 


6±4 

4 


0.29 ±0.03 

0.67 ± 0.13 




1.5 ±0.5 

5 


_2 

0.46 ± 0.05 




2.2 ±0.7 

6 


0.35 ±0.17 

0.3 ±0.5 





7 


0 ( fixed ) 

0 ( fixed ) 





2011 Outburst 

1 


_2 

1.09 ±0.08 


0.5 ±0.5 


9±3 

2 


_2 

0.4 ±0.3 




0.9 ±0.9 

3 


_2 

0.42 ±0.10 

5±3 

3.3 ± 1.3 


7.6 ± 1.0 

4 


_2 

0.60 ±0.18 

3.3 ±1.3 

1.6 ±0.8 


7.8 ± 1.2 

5 

0.50 ± 0.13 

0.54 ±0.12 

1.9 ±0.9 

4±4 


14 ( fixed ) 

10 ±5 

6 


0.23 ±0.04 

1.5 ±0.3 


0 ( fixed ) 


3±2 

7 


0.27 ±0.05 

1.0 ±0.2 




0.4 ± 0.2 

8 


0.27 ±0.05 

0 ( fixed ) 





9 


0.20 ±0.05 

0 ( fixed ) 




0 ( fixed ) 

10 


0 ( fixed ) 

0 ( fixed ) 
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Table 4 

RMS Amplitudes; r (%) 


Interval 

Lh2 

Lb 

Lh 

LIjY ^.^ow 

±hHz 

Li 

Lu 


1998 Outburst 

1 


14.4 ± 0.5 

17.8 ±0.7 




18.6 ± 0.7 

2 


19.0 ±0.3 

16.6 ±0.8 

16.8 ±1.6 



20.1 ± 1.1 

3 


19.9 ±0.4 

15.7 ± 0.9 

17.7 ±0.9 



20.9 ± 1.2 

4 


20.5 ±0.3 

15.5 ±0.6 

18.9 ±0.4 



20.0 ±0.7 

5 


20.2 ±0.7 

15.4 ± 1.3 

19.0 ±1.0 



19.2 ± 1.5 

6 


19.6 ±0.4 

17.5 ±0.7 

19.1 ±0.6 



19.1 ± 1.0 

7 


20.5 ±0.4 

15.7 ± 0.9 

18.3 ±0.7 



21.0 ± 1.1 

8 


19.3 ±0.7 

17.0 ± 1.3 

19.1 ±1.2 



19.0 ± 1.9 

9 


19.1 ± 0.5 

14.6 ± 1.5 

19 ±2 



16 ±2 

10 


18.6 ±0.6 

11.8 ± 1.7 

15 ± 3 



18 ±2 

11 


18.6 ±0.4 

10.7 ± 1.1 




24.2 ± 0.51 

12 


21.5 ±0.2 

8.2 ± 1.6 

19.0 ±1.8 



18 ±3 

13 


25.4 ± 0.3 


28.3 ±0.8 




2002 Outburst 

1 


17.7 ± 0.5 

13.7 ± 1.0 


8.4 ± 1.7 


8.0 ± 1.2 

2 


17.4 ± 0.2 

8.2 ± 1.6 

7.3 ±1.0 

11.8 ± 1.6 

4.0 ± 1.0 

8.7 ± 0.6 

3 

9.6 ± 1.3 

9.8 ± 1.3 

7.1 ± 1.1 

7.3 ± 0.9 

9.9 ± 1.1 


8.3 ±0.7 

4 

12.2 ± 1.0 

9.0 ± 1.1 

6.0 ± 0.5 

7.2 ± 0.4 

12.8 ±0.4 

2.5 ± 0.4 

7.6 ±0.6 

5 

6.8 ± 1.7 

14.2 ± 1.0 

9.3 ± 1.4 

7.8 ±1.3 

15.5 ±0.5 


9.7 ± 0.4 

6 


15.36 ±0.12 

15.3 ±0.3 


15.7 ± 0.7 


10.1 ± 0.7 

7 


15.05 ±0.13 

15.3 ±0.4 


15.1 ± 0.9 


10.6 ±0.8 

8 


15.17 ± 0.12 

15.6 ±0.4 


15.0 ±0.9 


10.0 ±0.9 

9 


16.08 ±0.18 

14.0 ±0.6 


17.4 ± 1.1 


10.4 ± 1.3 

10 


15.24 ±0.4 

14.1 ± 1.7 


15.6 ± 1.3 


10 ±2 

11 


15.7 ± 0.2 

14.0 ±0.7 


14.7 ± 1.2 


9.8 ±0.6 

12 


15.4 ± 0.3 

13.7 ± 1.0 


15.9 ±0.8 


8.9 ± 1.1 

13 


15.2 ±0.2 

15.9 ±0.6 


14.0 ± 1.7 


10.8 ± 1.2 

14 


14.7 ± 0.3 

14.5 ± 1.1 


15.7 ± 0.9 


9.2 ± 1.4 

15 


15.44 ± 0.19 

15.1 ± 0.6 


14.8 ± 1.1 


10.3 ±0.9 

16 


15.0 ±0.4 

17.5 ±0.7 


9.1 ± 1.5 


14.2 ± 1.2 

17 


15.4 ± 0.2 

15.7 ± 0.7 


14.5 ± 1.5 


9.8 ± 1.4 

18 


16.15 ±0.18 

15.3 ±0.5 


15.8 ± 1.6 


11.0 ± 1.8 

19 


15.7 ± 0.3 

15.3 ±0.6 


17.1 ± 1.3 


8.9 ± 1.6 

20 


15.6 ±0.2 

15.6 ±0.6 


15.2 ± 1.4 


9.2 ± 1.5 

21 


15.9 ±0.3 

15.0 ±0.7 


17.2 ± 1.3 


8.8 ± 1.9 

22 


15.7 ± 0.5 

14.3 ±0.9 




21 ± 2 

2005 Outburst 

1 


13.5 ±0.8 

18.4 ± 1.3 





2 


15.28 ±0.14 

14.5 ±0.4 


16.2 ±0.7 


7.4 ± 1.3 

3 


14.72 ±0.18 

13.0 ±0.6 


12.4 ± 1.0 


8.7 ± 0.7 

4 


16.47 ± 0.16 

11.9 ±0.5 


13.1 ± 0.6 


8.5 ±0.8 

5 


16.72 ±0.15 

11.6 ±0.6 


13.8 ±0.6 


8.8 ±0.5 

6 


16.6 ±0.2 

12.1 ± 0.7 


14.3 ±0.9 


5.9 ± 1.1 

7 

12.1 ± 0.9 

8.4 ± 1.2 

11.6 ±0.9 

5.0 ±0.7 

9.2 ± I. 4 I 


6.2 ±0.9 

8 


17.19 ±0.17 

12.9 ±0.5 


13.2 ± 1.0 


9.5 ± 1.0 

9 

8.7 ± 1.9 

13.2 ± 1.3 

14.1 ± 0.5 


12.8 ± 1.0 


10.6 ±0.8 

10 


15.8 ±0.2 

13.1 ± 0.6 


11.9 ± 1.6 


9.6 ± 1.4 

11 


15.9 ±0.2 

14.1 ± 0.8 


14.8 ±0.7 


8.6 ± 1.0 

2008 Outburst 

1 


15.0 ±0.2 

16.4 ± 0.5 


11.8 ± 1.2 


10.4 ± 1.7 

2 


15.5 ±0.3 

13.0 ±0.7 


19.1 ± 0.9 


9.1 ± 1.6 

3 


15.4 ± 0.2 

13.6 ±0.5 


15.8 ±0.7 


9.1 ± 0.8 

4 


15.3 ±0.5 

14.2 ±0.8 




12.3 ± 1.2 

5 


_2 

29.2 ± 1.0 




14.1 ± 1.6 

6 


11 ± 2 

14 ± 3 





7 


17.9 ±0.8 

23.5 ± 1.1 





2011 Outburst 

1 


_2 

20.2 ±0.9 


17 ±8 


11.5 ±0.8 

2 


_2 

14 ± 3 




13 ±4 

3 


_2 

23.7 ± 1.2 

7.2 ±1.4 

7.3 ± 1.1 


13.4 ±0.6 

4 


_2 

18.5 ± 1.7 

9±3 

9.2 ± 1.9 


13.4 ±0.7 

5 

6.8 ± 1.8 

14.3 ± 1.1 

7.8 ± 1.4 

5.3 ±1.2 


4.1 ± 0.8 

7.1 ± 0.9 

6 


17.0 ±0.5 

12.3 ± 1.2 


15.9 ± 1.7 


11 ± 2 

7 


16.5 ±0.6 

15.4 ± 1.2 




21.2 ± 1.8 

8 


16.4 ± 0.6 

21.4 ± 0.7 





9 


22.0 ±0.8 

18.7 ± 1.2 




28 ±3 

10 


33.8 ± 1.1 

29 ±4 
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Table 5 

Schechter function fit parameters. 


Interval a. 

^cut 

(Hz) 

(%) 


2008 

5 

-0.52 ±0.02 

3.33 ±0.10 

28.4 ±0.3 

2011 

1 

-0.450 ± 0.013 

5.63 ±0.09 

25.48 ±0.10 

2 

-0.79 ± 0.15 

5.6 ± 1.0 

15.4 ± 1.3 

3 

-0.59 ±0.02 

4.95 ±0.17 

26.75 ±0.15 

4 

-0.60 ± 0.02 

5.67 ± 0.16 

27.1 ± 0.3 


Over the following 14 days the flux and soft colors 
steadily decreased, while the hard color remained con¬ 
stant for 5 days before also starting to decay. The flux 
and soft color evolution transitioned into a faster decay 
on MJD 50929 (day 23). The flux decay briefly slowed 
down three days later. Finally the flux dropped below the 
detection limit on MJD 50934 (day 28), as the source re¬ 
turned to quiescence. During the entire decay SAX J1808 
remained in the extreme island state. 

4 . 2 . 2002 

In the 2002 outburst SAX J1808 was first detected with 
RXTE on MJD 52560. When pointed RXTE observa¬ 
tions commenced on MJD 52562 (day 6 ) the source was 
already close to the flux maximum. In interval I the 
power spectrum clearly reflects the island state, showing 
Lh, L/j, ThHz and In interval 2 both colors decrease as 
SAX J1808 appears to be transitioning to the lower-left 
banana. The power spectrum now shows a complex mor¬ 
phology, requiring at least six Lorentzian components. 
At ^50 Hz the Llf component appears, and we detect 
a 2.8cr feature at vq = 434 ± 12 Hz which we tentatively 
accept as the lower kHz QPO, as the upper kHz QPO 
is seen simultaneously at i/o,« = 597.9 ± 5.3 Hz. 

In the following observations a flux maximum of 
^85 mCrab is reached, and maintained for about a day. 
The lower kHz QPO is detected at higher significance 
with r'o,^ = 503.5 ±3.1 Hz and = 684.5 ± 4.9 Hz, 
and a second noise component appears near provid¬ 
ing further evidence that the source indeed transitioned 
to the lower-left banana. 

As the flux starts to decay on MJD 52564 (day 8 ), 
both colors increase and SAX J1808 transitions back into 
the island state. From interval 7 (day 9) the soft color 
starts to decrease once more, while the hard color re¬ 
mains roughly constant. The power spectrum shows very 
little change and remains consistent with an island state 
morphology. 

After reaching a flux minimum on MJD 52579 (day 23), 
SAX J1808 entered the prolonged outburst tail (Patruno 
et al.||2009), which continued for 24 days. 


4 . 3 . 2005 

Routine monitoring with RXTE revealed a new out- 
burst of SAX J1808 on MJD 53521 (day 1) (Markwardt 
et al. 2005). In the following days the flux steadily in¬ 
creased until a maximum of ~60 mCrab was reached ap¬ 
proximately five days after first detection. 


While the flux increased to its maximum (intervals 1- 
6 ) the hard color showed a steady decay, whereas the soft 
color remained roughly constant. During this time the 
power spectrum is that of the island state, showing Lf,, 
A/i, ThHz ^iud Ly^. 

The flux maximum of 60 mCrab is reached on 
MJD 53528 (day 6 , interval 7). In this interval both col¬ 
ors drop sufficiently to suggest a transition to the lower- 
left banana. The power spectrum confirms this transition 
with the appearance of ±62 and a narrow QPO at 58 Hz, 
which we identify as Tlf- Additionally a broad feature 
appears at ^'max = 428 Hz, which is difficult to identify 
as the low Q suggests it is ThHz, while the frequency sug¬ 
gests it is Li. The amplitude of 9% rms lies roughly in 
between LhHz and Li. We label it as ThHz, noting that 
it is not a ‘clean’ detection. 

After reaching the flux maximum, both colors in¬ 
creased as SAX J1808 transitioned to the island state. 
Over the following seven days the flux and soft color de¬ 
cayed, while the hard color increased. On MJD 53542 
(day 20) SAX J1808 was only marginally detected, mark¬ 
ing the end of the main o utburst and the start of a 40 day 
prolonged outburst tail (Patruno et al.|[2009). 


4 . 4 . 2008 

On MJD 54731 (day 1) SAX JI808 was again detected 
in outburst with RXTE. As in 2005 the rise to maximum 
flux was well sampled, and showed the same behavior: as 
the flux increased, the hard color decreased and the soft 
color remained constant. The power spectrum is an is¬ 
land state one and is well described with four Lorentzians 
(± 5 , Lfi^ ThHz and 7/^-) 

During its luminosity maximum (intervals 3-4, days 3- 
7) SAX J1808 reached a flux of 55 mCrab, the lowest 
flux maximum of all outbursts. The soft and hard colors 
again show a drop prior to the onset of the flux decay. 
Although this drop is as deep as seen in the other out¬ 
bursts, the power spectrum indicates the source does not 
make a state transition to the lower-left banana. A sharp 
feature does appear at iz = 283 Hz. Given = 520 Hz 
this is roughly where Li is expected to appear, but at 2a 
it is not significant. 

On MJD 54737 (day 6 ) both colors increase as the 
flux starts to decay. At this point (interval 5) the power 
spectrum shows an unexpected deviation from the regu¬ 
lar island state morphology as the high luminosit y flaring 
now dominates the power spectrum below 10 Hz (Bult & 
van der Klis 2014 1 . With a fractional rms amplitude of 

30%, the haring overwhelms the Lf, component that is 
normally seen in the same frequency range. We therefore 
cannot fit the Lf, component to the data. The flaring re¬ 
mains present for about three days, after which the power 
spectrum returns to the regular island state shape. 

On MJD 54748 (day 19) SAX J1808 enters the pro¬ 
longed outburst tail, which lasts for 30 days. 

4 . 5 . 2011 

The first pointed RXTE observation of SAX J1808’s 
seventh outburst was on MJD 55869 (day 2), roughly five 
days after t he initial detection with SWIET (Markwardt 
et al. The low frequency region of the power spec¬ 

trum was again dominated by high luminosity flaring. 
In the first two intervals the mid/high frequency domain 
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required three Lorentzians (L/i, LhHz and L„), which sug¬ 
gests that despite the low color values the source was still 
in the island state. 

On MJD 55871 (day 4, interval 3), a maximum flux of 
80 mCrab was reached, which correlated with a drop in 
amplitude of the flaring component. In the next interval 
the amplitude of the flaring increased again and the high 
frequency components suggested the source transitioned 
to the lower-left banana. The feature which is normally 
identified as the hump component has a central peak, 
but also broad wings and hence is not well described 
by a single Lorentzian profile, but instead requires two 
Lorentzians for a satisfactory fit: one to describe the 
sharp QPO and one for the overlapping noise component. 
We identify the broader noise component as Lh and the 
peaked QPO as Llf, but note that this identification is 
uncertain as the two components are not well resolved, 
and may be distorted by the high frequency wing of the 
high luminosity flaring component. 

In interval 5 (day 7) the high luminosity flaring dis¬ 
appears, while colors remain at low values. The Llf 
component is now clearly present in the power spec¬ 
trum and the marginal detection ( 2 . 8 cr) of the lower 
kHz QPO at vq i = 566.6 ± 9.4 Hz with a simultane¬ 
ous r'o,« = 736.1 ± 9.7 Hz confirms SAX J1808 entered 
the lower-left banana. 

In interval 6 (day 8 ) the colors increased as the source 
switched back to the island state. The power spectrum 
can be described with Lt, Lh, LhHz and L^. In the fol¬ 
lowing ten days SAX J1808 showed a regular decay with 
the soft color decreasing with flux while the hard color in¬ 
creased. During this time the power spectral components 
smoothly move to lower frequencies, before transitioning 
into the extreme island state in the last two intervals. 


4.6. The 410 Hz QPO 



the 3-60 keV energy range. We search all outbursts for 
this 410 Hz QPO using our 2-20 keV band, which gives a 
better signal to noise ratio. Due to the small amplitude 
of the QPO, we nearly always have to average multiple 
intervals to obtain a significant detection. 

We detect the 410 Hz QPO at > 3 ct significance in the 
outbursts of 2002, 2005 and 2008 (see Table |^, but not 
in the outbursts of 1998 and 2011, for which we obtain 
95% confidence upper limits for the amplitude of ~3%. 
The QPO amplitude varies between 2.1% and 3.6%, and 
shows no relation with ^’410 or iz„. The characteristic 
frequency varies smoothly between 409.6 and 416.6 Hz 
and is correlated with the upper kHz QPO frequency. 

We find the 410 Hz QPO is present in all observations 
in which has 320 < r'n< 401 Hz and the flux is higher 
than 40 mCrab (Figure]^. An exception is interval 14 
of the 2002 outburst (grey point in Figure zoom-in), 
for which we find a 95% confidence upper limit on r 4 io 
of 3.4%; comparable to the typical amplitude. Due to 
their lower count rates, the observations below 40 mCrab 
have higher upper limits on the amplitude (^4%) and are 
therefore also consistent with the Q PO be ing below the 
detection threshold (also see section |5.1.5|). 

The 410 Hz QPO is only seen in the island state, but 
is not related to a specific stage of the outburst evolu- 
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Figure 3. Overview of the 410 Hz QPO detections. We show flux 
vs. the upper kHz QPO centroid frequency for all intervals in grey. 
Intervals in which the 410 Hz QPO was detected are marked in 
yellow (2002), green (2005), and blue (2008). The vertical dashed 
line marks the spin frequency (t^spin = 401 Hz). The inset shows a 
zoom-in of the detection region. 

tion. In 2005 the QPO was detected during the rise to 
outburst, in 2008 at the maximum luminosity, while in 
2002 the QPO was seen during the flux decay. 

The proximity to the 401 Hz spin frequency suggests 
L 410 may be an upper sideband of the pulsations. If this 
is the case then a similar feature might be present at the 
lower sideband frequency of ~392 Hz. Using the same 
data selection for which we detected L 410 , we searched 
for the lower sideband, but no significant features were 
found, giving a 95% confidence upper limit on the frac¬ 
tional rms of 1.1% to 2.3%. T o improve sensitivity we 


use the shift-and-add method (Mendez et al. 1998) in 


which we shift the power spectra before averaging them 
such that the frequency of the predicted lower sideband 
feature, 2r'spin — 1 ^ 410 , is expected at 390 Hz. Again, a 
lower sideband is not observed. 

5. DISCUSSION 

It is common practice to study the variability proper¬ 
ties of X-ray binaries by considering the relation between 
the frequencies of several components of the power spec- 
trum. This approach is byed o n the succe ss of the WK 
( Wijn ands fc van der Klis||1999|) and PBK (|Psaltis et al.j 
999|) relations, which show frequency correlations over 


several orders of magnitude, linking the variability com¬ 
ponents of neutron star and black hole binaries. 

Applied to atoll sources the WK relation considers the 
correlation between the break frequency versus the Llf 
or Lh frequency, whereas the PBK relation considers the 
correlation between the low frequency QPOs and the sec¬ 
ond highest frequency feature in the power spectrum. For 
the higher luminosity lower-left banana state this would 
be L^, while for the lower luminosity extreme island state 
it is L^ow In the island state it is unclear which com¬ 
ponent should be conside red as neither L^ow or Li are 
observed (also see section 5.1. 2|). 


For atoll sources van Straaten et al. (2003) suggested 
a scheme of correlations tor the frequencies of several 
Lorentzian components plotted against the frequency of 
L„. This scheme of correlations encompasses the WK 
and PBK relations and is the same across different atoll 
sources. AMXPs also follow this atoll correlation scheme, 
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Table 6 

410 Hz QPO Fit Parameters 


Intervals 

1^410 

Q 

r 

Flap 

nb 

I'h 

I'u 


(Hz) 


(%) 


(Hz) 

(Hz) 

(Hz) 


2002 

6, 7, 8 

409.6 ±0.6 

80^29 

2.1 ± 0.3 

4.4 

3.07 ± 0.04 

18.30 ±0.10 

330 ±3 

10, 11 

414.7 ±0.6 

70l^® 

3.6 ±0.5 

4.3 

4.95 ±0.11 

25.3 ±0.3 

383 ±6 

12, 13, 15 

410.8 ± 1.1 

41+1’^ 

^-*^-12 

3.4 ±0.5 

3.6 

4.02 ±0.07 

22.1 ± 0.2 

349 ±7 

2005 

3 

416.6 ±0.8 

98t« 

3.2 ± 0.5 

4.1 

5.51 ± 0.09 

28.2 ±0.3 

390 ± 7 

2008 

1, 2, 3 

416.5 ± 0.9 


2.7 ± 0.5 

3.4 

4.55 ±0.08 

22.8 ±0.4 

399 ± 10 


Note. — Fit parameters for the 410 Hz feature, its significance (P/crp), and the simultaneously 
measured i^max of break and hump components and the upper kHz QPO. 


although they tend to be shifted in their by a constant 


factor of 1.1-1.6, which differs per source (|van Straaten 
et al.||2005 Linares et al.||2005 |. 


~K complication in considering frequency relations 
comes from the fact that the characteristic frequency 
of the Lorentzian function may not be the physical fre¬ 
quency of the variability mechanism. In this work we 
used iZmax as the characteristic frequency, which allows 
to treat QPO and noise components on the same footing, 
and tends to reduce t he scatter on frequency correlations 
(Belloni et al. 20021. However, models involving reso- 
nances or beat frequencies tend to be defined in terms 
of the centroid frequency vq = i^max/\/l + 1/(4Q^). For 
narrow QPOs (high Q) these two frequencies are almost 
the same, while for broad noise components (low Q) vq 
is significantly smaller than i^max- We will continue to 
consider the frequency relations in terms of i^max unless 
explicitly stated otherwise. 

Further difficulty in interpreting frequency relations 
comes from ambiguities in the identification of the vari¬ 
ability components. In particular when measurements 
deviate from a correlation, it could either indicate a break 
down of the relation, or a mislabeling of the considered 
component. It is therefore useful to distinguish between 
secure and more uncertain identification of the power 
spectral components. 


5.1. Frequency Correlations 

Figure 15 shows the frequency relations with respect to 
Vu for all outbursts of SAX J1808. Overall, the mea¬ 
sured frequencies follow the trends of the atoll correla¬ 
tion scheme. We quantify the relations by fitting a set of 
power-laws over the full Vu frequency range and present 
the fit parameters in Tablej^ We now discuss the notable 
features in this figure. 

5.1.1. Break and Hump 

The break and hump components are well correlated 
with Vu (see Figure®, but that correlation shows a break 
at i/u ~ 250 Hz. inis break could be due to a misiden- 
tification, however, as we are considering the three dom¬ 
inantly present power spectral components this is un¬ 
likely. 


The break in the correlation occurs at the transition 
between the island state and the extreme island state 
and appears to affect the break and hump components 
in exactly the same way. Indeed, if we consider the hump 
frequency as a function of the break frequency (i.e. the 
WK relation) we find the two components are tightly 
correlated (Figure and do not show state dependent 
behavior. This suggests that the break seen in Figure]^ 
is due to a change in that occurs when the source 
transitions between the island and extreme island states. 

When considering centroid frequencies instead of char¬ 
acteristic frequencies, the break remains present. Ad¬ 
ditionally, for centroid frequencies the WK relation no 
longer holds, with Lh drifting in frequency faster than 
Lb. 

If the upper kHz QPO indeed shows state dependent 
variations that do not occur in the other components, 
then considering frequency correlations against is not 
optimal. In Figure we therefore plot the frequencies of 
all components agains i//, (si milar to the i^and approach 
of van Straaten et al.||2002), which does not show this 
state dependent variation. 


5.1.2. Ltow and Lt 

The relation between L^ow and Li is a long standing 
problem. The two featu res were proposed to be due to 
the same component by |Psaltis et al.| (|1999D, who sug¬ 
gested that like the break and hump components, T^ow 
too becomes increasingly coherent for higher frequencies, 
eventually evolving into Plotting the frequency of Lh 
against that of Lfow or then reveals a tight correlation 
which can be extended over three orders of magnitude by 
linking black hole and neutron star data. 

The uncertainty of the PBK relation comes from the 
fact that the transition from L^ow to has never been 
observed. It has also been suggested tha t the two fea¬ 
tures are actually distinct components (van Stra aten 
et al.| 2003 2005). For SAX J1808 specihcally 


van 


Straaten et al. ] 2L)U5) suggested the two components were 


not the same because each requires a different frequency 
shift to be compatible with the atoll correlation scheme. 

Since most observations of SAX J1808 are in the is¬ 
land state, they do not have a L^ow or Li measurement. 
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Figure 4. Characteristic frequencies of the power spectral components as a function of i^u- Each point corresponds to a frequency 
measurement in Table The dashed lines show the fitted power-law relations (see Table and the vertical dotted line gives the frequency 
of the 401 Hz pulsations. For clarity we omitted points for which Qu was fixed at zero, as these points have large errors in Uu^ 
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Table 7 

Power-law Fit Parameters 


Component 

Data selection^ 

Normalization 

Index 

X^ldof 


Atoll Scheme Relations {ui vs i^u) 

break 

all 


(4.8 ±0.7) X 10“® 

2.29 ±0.02 

225IA7 

hump 

all 


(3.4 ±0.4) X lO-'^ 

1.86 ±0.02 

748/51 

hHz 


> 500 Hz 

(2 ± 5) X 10-6 

2.9 ±0.4 

6/7 

LF 

all 


(1.2 ± 1.0) X 10-6 

1.62 ±0.12 

6 .1/6 

b2 

all 


(5± 14) X 10-® 

3.3 ±0.5 

5.3/4 

law 

all 


(1±3) X 10-6 

3.4 ± 0.8 

6/9 

law Sz 1 

all 


(10 ±2) X lO-"' 

2.01 ± 0.04 

28/12 

LF & (410 - 401) 

all 


(1.2 ± 0.9) X 10-4 

1.97 ± 0.11 

58/11 

State Selected Relations {ui vs Uu) 

break b2 

EIS & all*^ 

(7± 2) X 10-6 

2.17 ± 0.05 

129/15 

break 

IS 


(2.1 ± 0.6) X 10-5 

2.05 ±0.05 

110/32 

hump 

IS 


(1.4 ±0.2) X 10-6 

1.64 ± 0.02 

519/35 

break 

EIS 

(0.4 ± 6.6) X 10-4® 

8.6 ± 2.6 

19/9 

hump 

EIS 

(0.2 ±3.5) X 10-4'4 

8.2 ± 2.5 

19/9 

WK relations vs uiy) 

break2 

all 


(10 ±6) X 10-2 

1.7 ± 0.2 

1.3/4 

LF 

all 


7.9 ± 1.0 

0.66 ±0.05 

6.8/4 

hump 

all 


6.6 ±0.1 

0.88 ±0.01 

186/54 

£ow 

all 


44 ±3 

0.60 ±0.09 

15/9 

u 

all 


208 ±2 

0.44 ± 0.01 

231/49 

LF & (410-401) 

all 


3.5 ± 0.7 

0.96 ±0.07 

83/9 

PBK relations {ui vs i^iowll) 

break 

all 


(1.2 ± 0.1) X 10-2 

1.15 ±0.01 

215/13 

hump 

all 


(1.12 ± 0.01) X 10-4 

1.04 ± 0.01 

148/12 

upper kHz QPO 

all 


31.2 ± 1.8 

0.50 ±0.01 

29/13 


Note. — Results of power-law fits to the frequency-frequency relations. 

IS — Island State, EIS — Extreme Island State 
^ i.e. EIS data for the break component and all data of the break2 component. 
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Figure 5. Characteristic frequencies of the power spectral components as a function of ui,. The green squares show the hump vs break 
frequency (WK relation). For further details see Figure]^ 
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Figure 6. Frequency measurements plotted against the L^ow 
(< 100 Hz) and Li (> 100 Hz) characteristic frequencies (PBK 
relation). Dashed lines show fitted power-law relations (see 
Table [^. 

For those observations that do have either component 
present we plot the measured frequencies against 
and vg in Figure The break and hump components 
and the upper kHz QPO all follow a correlation with 
viovi/i- Corresponding to the break below ~ 250 Hz 
observed in Figure the extreme island state Lu ap¬ 
pears to have a relatively high frequency, resulting in a 
flatter correlation for as compared to L^, and L^. The 
very similar correlations for Lt and Lh are in accordance 
with the idea that the T^ow and features are the same 
component, as assumed by PBK, and that it is that 
deviates at low frequency. 

Based on the PBK relation (Figure and Li 

could be the same, howeven when considering the fre¬ 
quencies against Vh (Figure™ it seems more likely that 
Lfow is somehow connected to LhHz- Whether or not the 
£ow component and the lower kHz QPO are the same or 
different is therefore not clear from the SAX J1808 data 
alone. 


5.1.3. The hHz Component 

The hHz component that appears in the power spectra 
of atoll sources is usually seen in a freque ncy range of 
100-200 Hz an d does not correlate with (|Ford & van 
der Klis 1998). Yet in SAX J1808 it appears th^ this 
hectohertz component does depend on (Figure]^. 

For Vu ^ 550 Hz LhHz appears in the same way as it 
does in other atoll sources, showing a large scatter be- 
tween 100 and 200 Hz, but rem aining constant with v,, 
( van Straaten et al.||2002 2005). Yet from — 550 Hz 


the frequency of the hectonertz component starts to sys¬ 
tematically increase with the upper kHz QPO frequency, 
in a way that is very similar to the trends seen in the 
other components. 

Systemat ic trends of z^hHz have p reviously been re¬ 
ported by Altamirano et al. (2008) who see a non¬ 
monotonic dependence of t'hHz on Vu- These variations 
appear similarly for > 550 Hz, yet are always within 
the scatter of i^hHz seen at lower i/„. In SAX J1808, 
however, the correlation extends well beyond the initial 
^'hHz scatter, suggesting the trend cannot be attributed 
to random variations in the frequencies. 

It should be noted that ThUz is not a dominant fea¬ 


ture in the power spectrum. Its properties may be influ¬ 
enced by the position and amplitude of the nearby kHz 
QPOs or potentially by an unresolved £ow component 
that could exist at similar frequencies. 

The correlation be tween AhUz and L„ was al ready ten¬ 
tatively observed by van Straaten et al. (2005), who also 
suggest it may be related to A^ow Indeed if the iow com¬ 
ponent follows the same frequency trend with as the 
break and hump components, then the island state ob¬ 
servations should show a blend of L^ow and LhHz- Then, 
as Vu increases, the iow component could become more 
prominent, causing the blended feature to follow the ap¬ 
parent correlation with iZu- 

If this interpretation is correct, it has a direct conse¬ 
quence for the relation between L^ow and L^. In one 
observation, interval 4 of 2002, the supposed £ow/hHz 
blend is seen simultaneously with a lower kHz QPO, 
which would imply that T^ow cannot be the same com¬ 
ponent as Lg. On the other hand, if Liow and Li are 
in fact the same component, like the PBK relation sug¬ 
gests, then the observed correlation between zzhHz and Uu 
might be intrinsic to the hectohertz component. Which 
of these interpretations is correct remains undecided. 

5.1.4. kHz QPOs 

So far, twin kHz QPOs had only b een reported once 


in SAX J1808 (Wijnands et al. 
spending to our interval 4 ol t 


2003|), from data corre- 
le 2902 outburst. This 


frequency pair is separated by roughly half the spin fre- 
quency and is consiste nt with ^ = |3r'o,^ — 2r'spin| 
( Wi.inands et al . |200 3) and a resonance at a 7:5 ratio 
( Kluzniak et al.|2l)U4 ). 


Unlike the t win kHz QPOs in most atoll sources (van 
der Klis|2006 ), for which the lower kHz QPO is the most 
prominent feature, in SAX J1808 it is the upper kHz that 
is most prominent. This ph enomenology is also seen in 
the AMXP XTE J1807-294 ( [Linares et al][2008| ) . 

In this work we present additional detections of the 
twin kHz QPOs in SAX J1808 in interval 2 of the 2002 
outburst and interval 5 of the 2011 outburst (Figure [^. 
The new detections of are only marginally significant 
(2.8 (t), however, their phenomenology is entirely con- 
sistent with the original detection by [Wijnands et al.| 
(20031, providing additional evidence that we indeed ob- 


serve the twin kHz QPO. All three observations are near 
the maximum luminosity of their respective outbursts, at 
a time when SAX J1808 has transitioned to the lower- 
left banana. The lower kHz QPO is very narrow and 
has a lower fractional rms amplitude than the upper kHz 
QPO. 

Considering the centroid frequencies of all three de¬ 
tections of the lower kHz QPOs we find that the twin 
kHz QPOs move over a frequency range of about 150 Hz 
(see Figure [8). This trend can be fitted with a constant 
frequency dmerence of i'o,u — i^o,e = (0.446 ± 0.009)r'spin, 
which is inconsistent with being exactly half the spin fre¬ 
quency, but, being slightly sm aller, may still be expl ained 
with a beat frequency model ( Lamb fc Miller[[2003 ). Al¬ 
ternatively the measurements may also be explained by 
a constant r'o,//t'o,u ratio of 1.36 ± 0.01, close to, but in¬ 
consistent with the proposed 7:5 resonance ratio. We do 
note that the obtained frequency trend is tentative, as it 
is based on a small number of measurements and may be 
subject to systematic uncertainties that are not captured 
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Figure 7. Twin kHz QPO’s in SAX J1808. The top panel shows 
the twin kHz QPOs in interval 2 of the 2002 outburst and the 
bottom panel those of interval 5 in the 2011 outburst. The pulse 
spike at 401 Hz was excluded from the fit. 
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Figure 8. Twin kHz QPO centroid frequency trend in 
SAX J1808. The dashed line shows the linear relation = 
Vu — 0.446^'spirn and the dotted line gives a proportionality at the 
constant ratio of 1.36. 
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5.1.5. The 410 Hz QPO 

Returning to the frequency relations shown in Figure]^ 
we now consider the ~410 Hz QPO. The characteristic 
frequency of L 410 correlates with the upper kHz QPO, 
although this relation is less s teep than the others. It 
was noted by Wijnands et al. (2003) that the difference 
between the 410 Hz QPO frequency and the 401 Hz spin 


frequency follows a quadratic trend with In Figure]^ 
we therefore subtract the spin frequency from the mea¬ 
sured L 410 frequency, revealing a frequency correlation 
for the 410 Hz QPO that is similar to that of the other 
components. This correlation suggests that the 410 Hz 
QPO is a sideband of the pulsations caused by the beat 
with a low frequency disk component. 

Because we consider rotation frequencies, a beat sig¬ 
nal may appear either as lower or upper sideband, de¬ 
pending on the relative sense of rotation. Since T 410 is 
observed as an upper sideband only, this suggests a beat 
with a quasi-periodic rotational low frequency disk phe¬ 
nomenon, L 410 - 401 , that is retrograde (counter-rotates) 
with respect to the neutron star spin. 

We searched for the associated low frequency L410-401 
component at 9 Hz in the power sp ectrum, both directly 
and using a shift-and-add method ( Mendez et al.|[l9M ), 
but did not detect it. Given that T 410 has a very small 
amplitude, it is possible that the original ~ 9 Hz sig¬ 
nal is below the detection limit in this frequency region, 
which is dominated by the prominent Lf, and Lh com¬ 
ponents. In the following we shall refer to the L 410-401 
component keeping in mind that this inferred ~ 9 Hz 
counter-rotational disk phenomenon was not, in fact, de¬ 
tected as a power spectral component, but seen only as 
410 Hz beat wit h the spin. 


It was noticed by van Straaten et al. (2005) that the 
T410-401 frequency is about hall Vh- 11 we fit the relation 
between these two characteristic frequencies we indeed 
obtain a proportional relation with a slope of 0.53 ±0.03. 
Using this relation we can predict the frequency at which 
±410 should appear, allowing for a shift-and-add search 
in the observ ation s where ±410 is not directly observed. 

In Section |4.6| we noted that L 410 is detected for 
320 < V < 401 Hz if the source flux is greater than 
^ 40 mCrab. We therefore first select all the low flux 
observations with in this frequency range for which 
±410 was not detected. A shift-and-add search revealed 
a marginal feature ( 2 . 8 cr) at Vspm-\-Vh/2, suggesting ±410 
might also present below ^ 40 mCrab. The observed fea¬ 
ture has a fractional rms of 1.4% for a fixed Q = 178. 
While this quality factor is higher than the typically ob¬ 
served Q ~ 80, it is consistent within typical errors. 

A shift-and-add search for all observations with i^o.u ± 
z/gpin and a search for all observations with i'q,u ± 320 Hz 
did not produce a detection (< Icr), giving 95% upper 
limits on the fractional rms of less than 1%. We note 
that this could either mean that the ±410 feature is not 
present in those observations, is too weak to be detected, 
or that the proposed relation with Vh is not correct. 

5.1.6. The LF QPO 

The narrow QPO at ^ 50 Hz, which we labeled ±lf, 


was identified as such by van Straaten et al. (2005) based 
on extrapolating the ±lf versus Lh relatioii round in low 
luminosit y bursters (IE 172 4-3045, GS 1826-24). It was 
shown by Altamirano et al. (2008) that the index of this 
power-law is consistent with 1 , implying the relation is 
consistent with a simple proportionality. Extending this 
dataset with the additional measurements from our work 
we find that a proportional function performs well and 
that the proportionality ratio is 0.65 ± 0.02. 

While the proportionality is consistent with 2/3, this 
small integer ratio of ±lf with Lh is only seen in 
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SAX J1808, IE 1724-3045, and GS 1826-24 (see Fig¬ 
ure ^ and only when considering characteri stic frequen¬ 
c ies Relat ivistic resonance models (Kluzniak & 


Abramowicz 2001), which use centroid frequencies, are 
therefore not applicable. 

While other atoll sources do not show a LF QPO that 
follows the same relation with Lh we see in SAX J1808, 
it has been suggested that some show a sub-harmonic 
iLi^F/ 2 , proportionality r atio of 0.33) rather than Llf 
( van Straaten et al.|2005|. How ever, many other sources 
have shown similar bJP Os (e.g. [Altamirano et ah 2005' 


20081 with frequencies drifting smoothly between Llf 


and Llf/ 2 j crossing the frequency range of the L 4 io_ 4 oi 
component (Figure!^. This suggests Llf, Alf/ 2 i and 
T410-401 might be tne same component. 



5 . 1 . 7 . LF and 4 IO Hz QPO 

If Llf and L 410-401 are the same feature then this 
raises the question why the LF component is seen directly 
in the power spectrum and L410-401 only as a sideband 
of the pulsations and why Llf and L 410 are never seen 
simultaneously. 

The 410 Hz QPO is only observed when the upper kHz 
QPO centroid frequency is below the spin frequency. We 
recently found that when the upper kHz QPO frequency 
moves above the spin freque ncy the pulse amplitude de - 


creases by a factor of ~ 2 (Bult & van der Klis 2015|). 
If this drop in pulse amplitude also applies to tne am¬ 
plitude of the beat signal, then the power of the 410 Hz 
QPO would drop by a factor of ^ 4, putting it below 
the detection threshold of those observations where the 
upper kHz QPO frequency is above the spin frequency. 

As the upper kHz QPO increases further, and presum¬ 
ably the inner edge of the disk moves closer to the neu¬ 
tron star surface, neither the LF QPO nor the 410 Hz 
feature are observed. Finally, when source transitions to 
the the lower-left banana, for the highest frequencies 
observed in SAX J1808 and thus when the accretion disk 
has probably advanced closest to the neutron star sur¬ 
face, the LF QPO appears. If the LF QPO originates 
from the inner region of the accretion disk, for instance 
through precession, then the increased disk temperature 
may be the reason that the LF QPO is observed directly, 
while the pulsations are still too weak to power a de¬ 
tectable beat signal. 

5 . 2 . Precession Models 

A popular interpretation of the frequency relations as¬ 
sociates QPO frequenc ies with relativistic orbital and 
precession frequencies (|Stella & Vietri||1998l |1999|). In 
this model the LF QPO frequency corresponds to the 
nodal precession frequency (or its second harmonic) of 
a test mass orbiting at the disk inner radius as caused 
by relativistic frame dragging, which in the weak field 
approximation is known as Lense-Thirring (LT) preces¬ 
sion. More recent versions of this model consider the LF 
QPO to be caused by a geometrically thick inner flow 
which, also due to frame dragging, precesses as a solid 
body (Fragile et al. ][2007l llngram et al.||2009|) , and pro¬ 
duces the cor rect frequencies m black note and neutron 
star systems (Ingram et al.||2009 Ingram & Done||20I0 ), 
while naturally explai niiig a raiige of QPO properties 
( Ingram fc Don^|20I2 ). 


Figure 9. Characteristic frequency of Llf versus for 
SAX J1808 and several other atoll sources. The solid line show 
a proportionality ratio of 0.66 and the dashed line shows a 
proportionality ratio of 0.33. 


For atoll sources specifically the relativistic precession 
models attribute Llf to LT precession. If L^p and 
L 410-401 are the same component, then the interpreta¬ 
tion that the 410 Hz QPO arises due to a retrograde 
beat with the neutron star spin, implies that the preces¬ 
sion too must be retrograde. However, L T precession is 
always prograde with the stellar rotation (Merloni et al.| 
1999 ] Markovic 2000), excluding frame dragging as the 


mam cause for the precession. 

For neutron stars, however, frame dragging is not the 
only cause of precession. The neutron star oblateness 
and its magnetic field may also apply torques to the in¬ 
ner accretion disk, which cause a retrogr ade (negative) 
contribution to the precession frequency (Lai||1999 Shi- 


rakawa & Lai 2002). In the magnetic precession model 


the stellar magnetic field induces a warp in a geometri¬ 
cally thin accretion disk, which in turn undergoes ret¬ 
rograde precessmn^jlo^yerjwith a power-law index of 
^ 4.6 (Shirakawa & Lai 2002), pure magnetic preces¬ 
sion predicts a much steeper relation with the upper kHz 
QPO than we observe. Classical precession (due to the 
oblateness), on the other hand, predicts a power-law in¬ 
dex of ~ 2.3, which is closer to what we observe. It would 
be interesting to investigate if in neutron star systems 
the combination of relativistic, classical and magnetic 
torque contributions can produce a retrograde var i ant of 
the LT thick flow precession seen by Fragile et al. (2007) 
in simulations of black hole disks, with all the attractive 
properties gained from it. 


6 . CONCLUSIONS 

We have presented a complete overview of the aperi¬ 
odic variability of SAX J1808 as observed with RXTE 
over the course of 14 years. We find that SAX J1808 is 
mainly observed as an atoll source in the island state, 
and sometimes transits to the lower-left banana or the 
extreme island state. The individual outbursts are very 
similar in terms of their light curve and color evolution 
and, aside from the 1-5 Hz high luminosity flaring ob¬ 
served in 2008 and 2011, show a similar evolution of the 
power spectra as well. 

Considering the power spectra of all outbursts we find 
that all characteristic frequencies are correlated with 
and that the correlations show a break when the source 
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state changes between the island and extreme island 
states around Vu — 250 Hz. This break disappears when 
considering the relations with Vb instead, suggesting that 
it is Vu that shows a state dependent relation with the 
other variability components. 

We find that at > 550 Hz the frequency of the hHz 
component in SAX J1808 shows an unexpected correla¬ 
tion with Vu, that roughly follows the frequency relations 
of the other components. We suggest that this correla¬ 
tion could be due to a blend with an unresolved low 
component or that it could be intrinsic to the hHz com¬ 
ponent itself. We also considered the relation between 
extreme island state £ow component and the lower kHz 
QPO, but found no definitive evidence in favor or against 
the two being the same. 

We also presented new measurements of the lower kHz 
QPO, which, for the first time, allows us to probe the fre¬ 
quency evolution of the twin kHz QPOs in this source. 
Our observations are consistent with a twin kHz QPO 
frequency separation of (0.446 ± 0.009)r'spin, which is in¬ 
consistent with half the spin frequency, but tentatively 
so, as the fit is based on a small number of measurements. 

Additionally, we presented the first detailed study of 
the 410 Hz QPO and found that it appears only when 
320 Hz < < Vspin- Subtracting the spin frequency 

from the measured 1 Z 410 we find that the resulting fre¬ 
quency matches the overall correlation trends, which 
strongly suggests that the 410 Hz QPO is caused by a 
retrograde beat against the 401 Hz spin frequency, even 
though the original ^ 9 Hz signal is never seen. 

Comparing the measured hump, LF and 410 Hz beat 
QPO with Llf detections in other neutron star systems 
we suggest that the LF and 410 Hz beat QPOs might 
be the same component. We suggest this QPO might 
be related to retrograde nodal precession caused by the 
(retrograde) classical and magnetic precession dominat¬ 
ing over the (prograde) LT precession. 

The authors acknowledge support from the Nether¬ 
lands Organisation for Scientific Research (NWO). 
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